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Abstract. We have obtained analytical exact solutions of the Dyer-Roeder equation in a cosmological model where creation of 
particles occurs at the expense of the gravitational field. We discuss the influences of inhomogeneities in the path of a light beam 
on the apparent diameter of astrophysical objects and consider both redshift independent and redshift dependent distributions 
of the inhomogeneities. 
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1. Introduction 

■ Cosmology has been , for a long time, a fertile ground for speculation. The choice between competing theories was very difficult 
I due to the small amount of reliable experimental data. 

Things have changed, however, in the last decades. The quantity of experimental results relating to the age of universe, its 
, expansion and its matter distribution, as well as gravitational lens occurrence statistics and related subjects, has grown to such a 
' extent that the room left today for speculative reasoning in Cosmology has been considerable reduced. 

Among the most interesting recent results, are the supernova lA type data, obtained at the end of the 1990s, which gave 
support to the hypothesis that our universe has an accelerated expansion CRiess. 1 998; Perlmuter. 1 998) . 

These observati ons lead to a revival of the cosmological constan t, as well as to new proposals for candidates able to generate 
a negative pressure llCaldwell et al.ll998HPeebles & Vilenkinll998 '). for example, quintessence. 

According to some of these hypothesis, the universe would have, beyond its usual baryonic matter content and dark matter, 
I also a negative pressure-generating con tent, a kind of dark energy that represents the vacuum contribution ( Peebles & Ratra 2002t 
i lCarroi200lUSahni & Starobinskvl2000l) . 

One of the attractive features of the hypothesis of particle production is that it therefore relates the large-scale properties of the 
universe to atomic phenomena ( Mac Crea 195 1 ). On the other hand, the introduction of this new component for the cosmological 
fluid gives not only an explanation for the cosmological acceleration, but also eliminates the age problem of the universe, which 
in the standard model is smaller than the one obtained for the age of the globular clusters. The estimate for the age of the 
universe depend upon the valu e of the Hubb le constant. If the value of Hq is near the upper limit obtained by Freedmann et al. 
(Ho = mkmlslMpc) jFreedman et al.ll997l) and consid ering the usual standard model (H - ^), we get some "relief for the age 
problem, although not a definitive solution (lOzeil 19991) . 

The model with particle production (OSC) provides also a reasonable fit with respect to kinetic tests, like luminosity di stance, 
angular diameter and the number counts of galaxies versus redshift relation jAlcaniz & Limall999tllTma & AlcaniJl999l) and in 
the radiation-dominated era (photon creation) the model can be compatible with present day isotropy and the spectral distribution 
of Cosmic Microwave Background Radiation (Lima 1997). 

The inclusion of A solves the ag e of the universe puzzle, but at the expense of creating a new one, the so-called cosmological 
constant problem dWeinherdfiW ). The conciliation between a very large value for this constant, predicted by quantum field 
theory, and a small one or zero, can be obtained if we consider the cosmological term time-dependent or quintessence models. In 
spite of, these models cannot explain why the dark energy density is comparable with the matter one. 

As an alternative model for the universe, we can introduce a cosmological particle production term, resulting in a scenario 
that can mimic the eff'ects generated by the inclusion of A. The physics involved is, nevertheless, quite different. 
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In this work we are going to study the exact solutions of the Dyer-Roeder equation, considering a homogeneous and 
isotropic universe where particle production occurs ait the expense of gravitational field energy {P rigogine 1989; Lima e t'al] 
Il996tlwaga & ArcurJl99'4llGariel & leDenmat ''l990: Zindhal 1993: Li ma & Alcanizll99"9tlAlcaniz & L ima 1999). 

In section 2 we outline the cosmological model and the relation between the source of particle production and the cosmo- 
logical term; in section 3 we examine the observational constraints for the model with respect to the age of the universe , the 
acceleration of its expansion, and make some considerations about primordial nucleosynthesis and structure formation, obtaining 
a narrower range for the values of the creation particle parameter; in section 4 we obtain the Dyer-Roeder equation for a universe 
with particle production and find exact solutions considering both a redshift independent and redshift dependent distribution of 
inhomogeneities . 



2. The model 

2.1. Field equations 

In this section we give a summary of the cosmology with particle production jPrigoginel 19891) , show how this model mimics the 
universe with a A term and find constraints in the creation parameter using the age of the universe and the deceleration parameter. 
The universe will be considered homogeneous and isotropic, described by the FRW line element 

ds^ = dt^ - R^{f){dp- + r^de^ + r^sin\0)dcf>^] . (1) 

The model with particle production is obtained from the field equations 

8;rGp = 3^+3^, (2) 
R K 

87r(P„, + P) ^ -2 T - -7 , (3) 

coupled to the balance equation for the particle number density 

h + 'i-n^'V. (4) 
R 

Here the dot means time derivative, p is the energy density, Pth is the thermodynamical pressure, and n and ^ are respectively 
the particle number density and the source of the particle production. Pc is the creation pressure, given by ( Calvao et al.- 19921) 

~P^_PlIlL^, (5) 

when the particle production is considered as adiabatic process. 

Combining equations (2), (3) and (4) with the state equation P,h = (v - l)p, it follows that 

RR.i^-p^)iR^..)^0, (6) 
2 2nH 

where v is a constant and k is the curvature. 

For a complete description of t he cosm ological environment we must have an explicit expression for the source of particle 
production. Following Lima et al. jLima et al. 1996), a physically reasonable expression for the particle creation rate is *P = 
3nfiH. The creation parameter /3 is generally a function of the cosmic era, or equivale ntly of the paranieter v . We can assume 
two different parameter creations, /3y and p,„, for radiation and matter eras, respectively caniz & Lima '199^. However, since 
we study a distance treatment in the actual matter dominated phase (v = 0), we supposed /3m = /3, with /3 constant. Beyond this, /3 
must be positive, otherwise the second law of thermodynamics iS"-a > 0) is violated. 

Substituting the particle production rate above in equation (6) and integrating , we obtain 

R = /;o(f )^ , (7) 

ro 

where the subscript refers to the present time. 

The OSC can mimic a universe with a A term. To clarify this relation we take into account the entropy production in both 
cases. For OSC we have 

S.„ = na- + crip = -— , 

while for a universe with a A term 

= ^ 
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T is the temperature, yu = - Tcr is the chemical potential and cr is the entropy per particle. 

The two scenarios can explain the entropy co ntent of the univ erse, which is poorly explained using only Einstein's field 
equations, which are invariant under time reversal (Prigogine ^l989l) . 

Considering the process to be adiabatic (& - 0) and equalizing the specific entropy due to the A term and due to particle 
production we obtain 

— oc ^ 

SttG n 

For ij/ - 3npH this implies A oc r^, for a pressurelles universe. The two scenarios have identical time dependence for the 
scale factor, and consequently for the Hubble function, too. 

2.2. Observational constraints 
2.2.1 . Age of the universe 

Using the scale factor (7) the Hubble function is given by 

^= ^[(l-yS)f]-'. (8) 

The results from the Hubble Space Telescope Key Project (lFreedmanl2000l) indicate. 
Km 

//o = 72±8- , (9) 

Mpc ■ s 

for the Hubble constant. 

T he High Z Supernovae estimates for the age of the universe is 14.2 + 1.7 Gyrs (lRiesJl998l) and 14.9 ±1.4 Gyrs jPerlmuteJ 
Il998 >). Carreta et al.(2000) using results from Hipparcos, RRLyrae and Chepheids to re-calibrate the globular cluster distance 
scale, found 12.9 ± 2.9 Gyrs. Butcher (1987) used the abundance of the Th^^ to estimate the age of the stars and consequently a 
lower limit for the age of the universe. The interest in Th abundance is related to its half life, about 14.05 Gyr, an interesting time 
from a cosmological point of view. Using the relation between the abundance of Th and Eu, Westin et al. (2000) estimate for the 
universe an age of 15.0 Gyr, while Johnson and Bolte (2000) found 1 1.4 Gyr Recently Krauss and Chaboyer (2000) estimates a 
lower limit on the age of the universe, 1 1 Gyr with 95% confidence. 

These values differs substantially from each other, depending on the technique used; we update the interval for the value of 
the /3 parameter found by Lima et al. (Lima et al. 199a). 

Using (8), the upper limit for //q from Freedman et al. and the estimated age for the oldest objects in our galaxy {l6.3Gyr), we 
obtain the value 0.49 for the upper limit for the creation parameter To obtain the lower limit for /?, we use the same expression 
(8), but with the lower limit for Ho, again from Freedman, and the smallest age for the local oldest stars from Carreta et al. (2000), 
giving 0.02. 

So, the validity interval for the creation parameter is 
0.02 </3< 0.49 . 



2.2.2. Acceleration of tine universe 

According to data from type la supernovae observations, the universe is accelerating. Generally, this acceleration is thought to 
be due to repulsive gravity, which can be introduced via a negative pressure, a cosmological term or a quintessence component 
in the matter content of the universe. 

An alternative view would be to consider particle creation, which naturally redefine the energy-momentum tensor and could 
account for this increasing expansion velocity ( Freaza et al. 2002; Alcaniz & Lima 1999; Lima et al. 1996; Zindhal 2002J. 

Taking into account (6), we can write the deceleration parameter in terms of the particle creation source, namely 

1 A-TTmCV 

G = (10) 

where m is the rest mass of the produced particles, and we consider a null value for the curvature. We obtain an accelerated 
universe, expanding with the scale factor (7), if/? > ^. Taking into account the age of the oldest objects in our galaxy, the 
estimate for the Hubble constant and the indications that we live in an accelerating universe, we obtain that 

0.33 < yS < 0.49 . 
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2.3. Primordial nucleosyntliesis 

In spite of the cosmological consequences, t he curr ent matter creation rate , ~ 10"'^ nucleons cm^^yr^^ , is nearly the same 
rate predicted by the stea dy-state Universe jHovld i 1993) taking into account k - 0. This matter creation rate is presently far 
below detectable limits Iai caniz & Li ma 1999). Consequently, this indicates, probably, that the early picture of the OSC differs 
slightly from the standard model with respect to the primordial nucleosynthesis. However a quantitative analysis with a multifluid 
description is necessary to more accurately determine the primordial nucleosynthesis in OSC. 

Considering primordial nucleosynthesis in its standard form with 3 neutrinos flavours, the main parameter of the theory is the 
baryon number density quoted as its ratio to the number density of photons, namely t] - j-. Note that, in a given era of evolution 
of the universe, the model with cosmological particle production affects the number density of photons in the same way that the 
number density of the baryons. So, in this model rj does not suffer sensitive variations in relation to the standard model. 



2.4. Growth of structures 

In the OSC model the particles are created from the gravitational field energy. In this case, a negative extra pressure appears in the 
energy momentum tensor This extra pressure is responsible for cosmic acceleration and has the same influence as the negative 
pressure that appears in models with a cosmological term. So, in the same way that the effects of a cosmological term on linear 
perturbation theory over short distances can be studied using the expression 

6 + 2H6 - AnGpd = , 

we consider that this equation is also applicable to the OSC model, where 6 is the density contrast. Substituting the scale factor 
(7) in the equation for the density contrast evolution we obtain the following modes 

6- oc t f-'fi-') 



For p = we obtain the usual growing and decaying modes for the dust Friedmann model, 5+ oc ft, 5_ oc r\ respectively. Note 
that the 6+ for the OSC model grows faster than the growing mode without particle production (J3 - 0), including for accelerated 
expansion of the universe. The few indications in this section are, from our point of view, sufficient to motivate a more detailed 
study of the scalar density perturbations to furnish a more precise picture of the growth of structures in the OSC model. 



3. The Dyer-Roeder equation 

If there is matter distribution along a cone defined by the bundle of light rays connecting the source to observer, the angular 
diameter distance of that source from the observer is smaller than that which would occur if the source was seen through an 
empty cone. This effect is known as Ricci focusing. Our intention in this section is to find and solve the mathematical expression 
of this effect, the Dyer-Roeder equation, when the particle number in the universe is not conserved. 
For a conformaUy flat metric the Ricci focalization equation is 

Va + F Va = , (ii) 

where the dot means derivative with respect to the affine parameter, A is the beam cross sectional area and - ^ is the vector 
tangent to the photon trajectory ( lDemianskil2000l) . Considering an energy momentum tensor of the form 

T^V^{P + P)u^Uv + Pgtiv , (12) 

where the tot al pressure P inclu des contributions from the thermodynamical and creation pressures, we obtain, after straightfor- 
ward algebra JGiovi et alJ200(]t) 



Va + 47r^^(i + z)2 Va = . 



^0 



(13) 



Taking into account equations (5) and (7), we obtain, after a simple but a lengthy manipulation jGiovi et alJl2000l) . the Dyer- 
Roeder equation in the cosmological background with particle production 

(1 + ^)'tT + ^^(1 + ^) T + -/?)[1 = 0, (14) 

flz^ 2 dz 2 

where we have defined the dimensionless distance r - DHq, D as the angular diameter distance. The paramet er a is related to 
the R icci focusing and conceptually is the fraction of homogeneously distributed matter inside the light beam jDver & RoedeJ 
Il974l) . 

Let us consider two interesting particular cases. 
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3.1. Constanta 



We recall that, for a = 0, our model represents a universe where all matter is clustered, while, for a = 1, we have a universe with 
its matter content homogeneously distributed. 

Equation (14) can be analytically solved for any value of a in the interval from zero to one. We obtain 



riz) = (1 + zr^'\CiJ^ [A] + C2Y^ [A]} , 



where 



A = ( 



(15) 



(16) 



J and Y are the first and second Bessel functions, respectively. 

The profile for r{z) when the universe is non accelerated (jS < j) is shown in Fig. 1. 




0.2 0.4 0.6 0.8 1 1.2 1.4 , 1.6 1.8 2 2.2 2.4 2.6 2. 



Fig. 1. Evolution of Dyer-Roeder distance for j6 = 0.1, and three different values of a. 
For an accelerated universe (jS > i) the profile is the one seen in Fig. 2. 
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Fig. 2. Evolution of Dyer-Roeder distance for j6 = 0.4, and the same three different values of a. 

To obtain the limiting case when there is no particle production we must put y8 = in the original differential equation (14). 



3.2. Variable a 

In this section we solve the Dyer-Roeder equation for the case where a is no longer a constant, but instead a particular function 
of the redshift. We adopt for aiz) an expression identical to the one used by Linder (1998), namely 



a(z) = p{l+zf, 



(17) 



where p and q are constants. It is reasonable to expect the clumping to decrease with increasing redshift, implying ^ > 0. Some 
care is needing in choosing the parameters p and q is due to the range of a, that must stay in the interval < or < 1. Naturally, 
this will depend on the range of redshift studied. 
Integration of equation (14) now results in 



riz) = .ClJA(B) + .C2YA(B), 



(18) 



6 

where 
A 
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2 iS-Zp + q) 



and 



B = 



The profile of r{z) in this case, is given by Fig. 3, for specific values of q and p, when the universe is not accelerated (J3 




0.4 z 0-6 



Fig. 3. Evolution of the Dyer-Roeder distance for ji - 0.25, and two different values of q and p. 
For an accelerated universe (fi> the profile is illustrated in Fig. 4. 



r(z) 




Fig. 4. Evolution of Dyer-Roeder distance for jS = 0.8, and same two different values of q and p. 
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4. Conclusions and final remarks 

The cur rent matter creation rate, *Po ~ 10"'^ nucleons cm^^yr^^ , is nearly the same rate predicted by the steady-state universe 
reg ardless of the value of the curvature parameter. In spite of the matter creation rate presently being far below 
detectable limits ( Alcaniz & Limall999l) . the physical consequences for the evolution of the universe are potentially measurable. 
We have updated (and narrowed) the range of allowed values for the creation parameter 

For the Dyer-Roeder equation, we have solved it in a background universe with particle production, expanding considerably 
th e manifold of solutions. As far as we know, this equation has previously been solved only for some special constant values of 
ryllrTiovietaU2nnn). 

In Fig. 1 we have plotted r{z) for a non accelerated universe, for three different values of constant a. We note that increasing 
the inhomogeneities in the path of the light beam (smaller a) leads to greater values for the Dyer-Roeder distance at the same 
redshift. Remembering the definition of r(z), this implies smaller values for the viewing angle as a diminishes. 

We can interpret these results saying that the objects will appear at apparent distances greater than the ones we would observe 
in a totally homogeneous universe (a = 1). 

In Fig. 2 we repeat the same graph, but now in an accelerated universe. We observe the same qualitative behavior, with little 
enhancing of the effect. The origin of this effect is in the inhomogeneities, while the acceleration, although having some influence, 
is not the main agent in this case. 

In Figs. (3) and (4) we show r(z) in the case of variable a, given by (19), for non accelerated (fi - 0.25) and accelerated 
(jS = 0.8) universes, for different values of q and p. In both cases, for the redshift interval studied, the range of the smoothness 
function a{z) lies in the interval < a(z) < 1 . 

The next steps for future work are to analyze the lensing probability occurrence for OSC. 
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